Stimulated Neutrino Transformation in Supernovae 
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Large amplitude oscillations between states of a quantum system can be stimulated by sinusoidal ex- 
ternal potentials with frequencies that are equal to the energy level splitting of the states or a fraction 
thereof. We investigate this effect for neutrino oscillations both analytically and numerically finding a 
simple expression for amplitude and wavelength of the transitions as a function of the density, the ampli- 
tude and wavenumber of the fluctuation, and the matter mixing matrix elements. We apply our findings 
to the supernova environment and find that it is possible to predict stimulated transitions that occur due to 
sinusoidal perturbations in the density. 

PACS numbers: 14.60.Pq,97.60.Bw 



INTRODUCTION 



Neutrino flavor transformation is a complex phe- 
nomenon: the general coherent neutrino state is a linear 
combination of several states with different energy and the 
oscillations between these states are affected by the pres- 
ence of matter and may exhibit many-body effects, e.g. 
OJj, if the neutrino density is sufficiently high. Yet de- 
spite this complexity, neutrino flavor transformation is a 
quantum phenomenon so that one should also expect an 
additional type of transformation: stimulated transitions 
driven by perturbing potentials with appropriate frequen- 
cies. This possibility of driving transitions between the 
neutrino states was recognized many years ago, e.g. 
when periodic potentials were studied for their relevance 
for neutrinos passing through the Earth and the Sun. It 
was found, e.g. |5|], that large amplitude transitions - called 
parametric resonances - occur between the neutrino states 
when the frequency of the perturbation, k+, matches an in- 
teger fraction of the frequency associated with one of the 
energy splittings, i.e. nk+ = 6k a b. This result is intrigu- 
ing because any density profile can be decomposed into a 
Fourier series, even non-periodic potentials e.g. turbulence 
10,01, and some terms in the series may have wavelengths 
which are capable of driving parametric resonance transi- 
tions. But even if such a mode is possessed by the profile, 
in order to observe the transition the neutrino must traverse 
a distance 1njK n , where 1/k„ is the reduced wavelength of 
the stimulated transition at the parametric resonance. This 
wavelength is a function of both the amplitude of the den- 
sity perturbation, i.e. the Fourier coefficient, the perturba- 
tion wavelength 1/&* and the integer n. But one needs rela- 
tively high densities and/or large amplitude fluctuations for 
an appreciable effect to occur over a reasonable distance 
and such conditions are not found in the Earth or in the 
Sun. 

The combination of large amplitude (~ 10%) density 
fluctuations and high densities can be found in supernovae 
and supernova density fluctuations and were found numeri- 
cally to have a large impact on neutrino flavor evolution [8|- 



fToh - With this motivation, in this letter we revisit the phe- 
nomenon of stimulated transition for the supernova envi- 
ronment. We first derive expressions for the reduced wave- 
length and the amplitude of stimulated transitions both on 
and off the parametric resonance. Next, to illustrate our 
findings, we examine the case of a specific supernova den- 
sity profile which includes a forward shock to which we 
add a perturbing sinusoidal fluctuation. Using this sce- 
nario, we predict the places where the transitions will oc- 
cur and then verify our predictions with a numerical three 
neutrino oscillation calculation. Our example demonstrates 
that it is possible to identify fluctuation wavelengths that 
will induce large amplitude transitions in the supernova en- 
vironment. 



STIMULATED NEUTRINO TRANSFORMATION 

We are interested in the probability that some initial neu- 
trino state \v(x)) at x is later detected as the state |v(V)) at 
x' . These probabilities are computed from the 5 -matrix 
which relates the initial and final states via \v(x')) = 
S (Y, x) \v(x)). The probabilities we calculate depend upon 
the basis and to avoid the intrinsic oscillatory behavior in 
the flavor basis we prefer to use the instantaneous eigen- 
states of the Hamiltonian H, known as the matter basis. 
The two bases are related by the unitary transformation 
matrix U. For three flavors/mass states the matrix U is 
parametrized by three mixing angles, 612, #13 and 623, a CP 



phase and two Majoranna phases. The differential equation 
for the neutrino S -matrix jfTTl[T2h is simply 



ax 

and in the flavor basis the Hamiltonian is 



//(/) = U K™UZ + 



(1) 



(2) 



where is a diagonal matrix in the vacuum, Uq the 
vaccum mixing matrix, and V^\x) some 'potential' that 
we allow to be position dependent. Our goal is to solve 
Eq. (Q]) for S given some potential V(x). The problem we 
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have in mind is the case where V(x) possesses some sort 
of 'smooth' component we denote by V and a perturbation 
5V. Thus the flavor basis Hamiltonian is split into three 
terms, H (f) = U U f + V if) +5V (f \ which we group as 
ffC/J - fj<j) + SV^K The matter states are defined as those 
which diagonalize i.e. if is the diagonal matrix of 
eigenvalues of H then U is = U K (m) Similarly, 
the unperturbed matter states are defined as those which 
diagonalize H (f) i.e. H (f) = U K m t/ 1 " where K (/h) is the 
diagonal matrix of unperturbed eigenvalues, ki,k,2, . . ., of 
H. In this basis 



H' 



(in) _ iH/ft) 



dx 



(3) 



We now write the S -matrix for the matter basis as the prod- 
uct S (m) = S A where S is defined to be the solution of 



dS_ 

dx 



dU 



dx 



(4) 



If we know the solution to the unperturbed problem, S , we 
can solve for the effect of the perturbation by finding the 
solution to the differential equation for A: 



dx 



(5) 



In general the term wSV^'U which appears in this equa- 
tion possesses both diagonal and off-diagonal elements. 
The diagonal elements are easily removed by writing the 
matrix A as A = W B where W = exp(-iH) and S a diago- 
nal matrix S = diag{^\,^ ■ • ■)• Substitution into equation 
d5j gives a differential equation for B 



i— = W 

dx 



dx 



WB 



(6) 



and S is chosen so that dE/dx removes the diagonal ele- 
ments of S^WSV^'U S . Once S has been found, deter- 
mining transition probabilities involves solving dB/dx. 



CONSTANT POTENTIALS WITH SINUSOIDAL 
PERTURBATIONS 

We now consider the specific case of a matter potential, 
so that the only non-zero entry of is V ee and V ee (x) 
is V ee (x) = vIGf h e (x) where h e is the 'smooth' electron 
density. With this form for the potential the elements of 
f/tyWf/ = V m are = U*U ej V ee . Next we Taylor 

expand V ee (x) as V ee (x) = V* + . . . and retain only the first 
term of the expansion. For constant V ee (x) = V+ the S- 
matrix describing the unperturbed problem, S , is diagonal 



because the eigenvalues of H and thus the mixing matrix U 
are fixed so S is simply S = exp (-iK^ x). The perturba- 
tion we consider is the case of a single sinusoidal fluctua- 
tion of wavenumber amplitude C* and phase shift x* 
i.e. 8V ee (x) = C* V* sin [k+(x + x*)]. With this perturbing 
potential we solve for the quantities 

6 = C * V * ^ [cos (k+x+) - cos (fc* [x + x*])] . (7) 
k* 

Everything we require to start solving the equation for 
dB/dx has now been defined. For just two flavors the dif- 
ferential equation for the matrix B is 

dB 

i— = C*V* sin(fc* [x + x*]) 
dx 







(8) 



e2 







B 



where 5k\2 = k\-k2 and similarly d^n = g i - £2- Our first 
step in solving this equation is to introduce the quantities 
Z12 and K n , defined to be 



Z12 



(|t4l| 2 -|^2| 2 ), 



{-I) : Jn{Zl2)U el U e 2 



Xe 



Z12 



(9) 



(10) 



respectively where /„ is the Bessel J function, and our sec- 
ond is to use the Jacobi-Anger expansion for exp (idgu) 



exp (16^12) = exp[izi 2 cos(fe*x*)] 



^ {-i) n J n (Z12) exp [1 nk+(x + xj] 



(11) 



With this expansion and two definitions we find equation 
([8]) becomes 



dB 

dx 



= iY( ° 

/ l y K * e -i\SkYi+nK\ 







B (12) 



which is a greatly simplified equation. We now adopt the 
Rotating Wave Approximation and drop all terms in the 
sum except the one closest to parametric resonance i.e. 
where 8k\2+n w 0. With this approximation we find that 
equation ( TT2T ) falls into a class of matrix Schrodinger equa- 
tions with known solutions. Introducing 2k„ = Skn +nk+ 
and q\ = k\ + k„, our final result is that 



B 



( e ik n x [ cos (q nX ) _ ^ sin(^r„x)] 
e~ lk " x ^sin(q„x) 



-ik n x 



_ e ,k„x ^ s j n ( ?nX ) 

[cos(^r„x) + 1^ sin(^r„x)] t 



(13) 
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FIG. 1: Using the base density profile described in the text, we 
plot multiples of the wavelength associated with the mass split- 
tings, n/dk\2 (upper panel) and n/dkii (lower panel) as a func- 
tion of distance. In each panel, the lowest solid curve repre- 
sents n=l, the middle is n=2 and highest is n=3. The wave- 
length of a sinusoidal perturbing potential with reduced wave- 
length /t* = l/ki, = 10 km is plotted as the dashed line. When 
any of the solid lines crosses the dashed line, the parametric res- 
onance condition is fulfilled and thus there exists possibility of a 
stimulated transition. 



Because both S and W are diagonal matrices, from this 
equation we can read off the transition probability between 
the matter states 1 and 2 
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P i2 = sin 2 (q n x). 



(14) 



The range of k+ such that the amplitude of the oscillation is 
greater than 1 /2 is = 4|/c„|/« and we see that this width 
is proportional to C™ when the parameter z.n is small. In 
summary then, we expect to see in Pu(k+) a series of peaks 
at the undertones of [Sk^l with widths that decrease geo- 
metrically. This same behavior has been found for irradi- 
ated two-level molecules with permanent dipoles IU3I1 . 



APPLICATION TO THE SUPERNOVA ENVIRONMENT 

We demonstrate the phenomenon of stimulated transi- 
tion in supernova-like conditions by constructing a density 
vs. radius relationship from the parametrized form from 
Fogli et al. IToh . The shock position was set at t = 4 s 
and we superpose on the profile a single sinusoid perturba- 
tion with a reduced wavelength /t* = = 10 km and 
amplitude C* =0.1. Through this profile we send a neu- 
trino with the following properties, d>m 2 2 = 3 x 10~ 3 eV 2 , 



Using the base density profile and the parametric reso- 
nance condition, Skjj + nk+ = 0, where the subscripts i 
and j denotes matter states 1,2 or 3, we predict the loca- 
tions where significant stimulated transitions between the 
states may occur. These are the places where the n/Skij 
curves in Fig. ((TJ, intersect the horizontal dashed line cor- 
responding to X*. We see that there are seven locations 
(c,d,f,g,i,j,m) in the 1,2 channel (top panel) and seven loca- 
tions (a,b,e,h,k,l,n) in the 2,3 channel (bottom panel). 



5m\ 3 = 8 x 10" 5 eV 2 , <9 12 = 33°, 
and with an energy of 20 MeV. 
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9°, 023 = 45°,£ = 



The formalism developed in the previous section as- 
sumed a constant density profile, but the supernova-like 
density profile we have chosen is not constant. Thus we 
only expect to see stimulated transitions at some subset 
of the fourteen locations identified in Fig. [TJ To deter- 
mine this subset, can use a comparison of the density scale 
height r p = \p/(dp/dr)\ with the reduced parametric reso- 
nance wavelength l/q n . When r p » l/q n , then the con- 
stant density approximation is reasonable, and we should 
expect to see stimulated transitions. In Fig. [2] where the 
top (bottom) panel again corresponds to the 1,2 (2,3) chan- 
nel, it can be seen that the stimulated transition wavelength 
becomes long at the parametric resonance locations. In the 
1,2 channel, only for points f,g is r p » l/q n , although 
r p is slightly larger than l/q n at m. Thus, we expect a 
significant stimulated transition at f,g, and a more mod- 
erate one at m. In the 2,3 channel, this figure suggests that 
the most significant transition will be at location k with 
more moderate transitions at h and 1. In Fig. [3] we check 
these results by evolving the neutrinos through the base 
+ perturbing profile with an exact numerical three neu- 
trino calculation. The base density profile possesses three 
Mikheev-Smirnov-Wolfenstein (MSW) [Qj, [H high (H) 
density resonances. Two of these, labelled X and Y on the 
figure, are semiadiabatic while the third, Z, is at the po- 
sition of the shock and is diabatic. The evolution of the 
neutrino through the three MSW resonances is accounted 
for with the unperturbed matrix S . Note that we have plot- 
ted transition probabilities between matter eigenstates not 
mass eigenstates, so the relevant MSW H transition is in 
the 2,3 channel. At the expected stimulated transition lo- 
cations, we indeed see strong transitions at f,g between 
matter eigenstates 1 and 2, and a somewhat smaller tran- 
sition in the same channel at m. In the 2,3 channel we 
see a strong transition at k and smaller transitions at h and 
1. We see also a few additional features at the parametric 
resonance locations caused by simultaneous occurrence of 
stimulated transitions between multiple pairs of states. For 
example, around points f,g and around m both 1 /6ki2 ~ /t* 
and 3/5^23 ~ This dual satisfaction of the paramet- 
ric resonance conditions implies 4/5k\3 ~ A* which means 
we can have stimulated transitions between all three states 
simultaneously. 
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FIG. 2: Comparison of the reduced wavelength 1 /q„ for the n=l, 
n=2 and n=3 modes (red, blue and purple lines) with the density 
scale height, r p for the base density profile (black dashed line). 
At the parametric resonances identified in Fig. Q] the transition 
wavelength \/q„ becomes long. When r p is larger than l/q„, we 
expect stimulated transitions. 



CONCLUSIONS 

We have investigated the effect of sinusoidal density 
fluctuations upon neutrino propagation and derived an ana- 
lytic solution for the matter basis transition probabilities in 
the case of constant density. Large transitions between the 
states can occur and we have determined their wavelengths 
and amplitudes. We presented an example calculation with 
a supernova-like base density profile to which we added a 
sinusoidal perturbation. When compared with a numerical 
three neutrino flavor calculation, we found we were able 
to successfully predict the locations of the stimulated tran- 
sitions, despite the presence of multiple MSW resonances 
in the unperturbed profile. This method holds promise as 
a tool for analyzing consequences of density perturbations 
on supernova neutrino flavor transformation. 
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02ER41216 (GCM+KMP) and by a NC State University 
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FIG. 3: Transition probabilities between matter eigenstates, P\2, 
Pi3 and P23, as calculated numerically. Large stimulated transi- 
tions occur at f and g (top panel) and at k (bottom panel). Smaller 
amplitude stimulated transitions occur at m, h, and 1. The tran- 
sitions labeled X, Y and Z are base profileMSW H resonances. 



[1] H. Duan and J. P. Kneller, J. Phys. G 36 113201 (2009) 
[2] Schafer, A., & Koonin, S. E., Phys. Lett. B 185 417 (1987) 
[3] Krastev, P. I., & Smirnov, A. Y., Phys. Lett. B 226 341 
(1989) 

Koike, M., Ota, T, Saito, M., & Sato, J., Phys. Lett. B 675 
69 (2009) 

Akhmedov, E. K., Nuclear Physics B 538 25 (1999) 
A. Friedland & A. Gruzinov, astro-ph/0607244. 
[7] J. P. Kneller and C. Volpe, Phys. Rev. D 82 123004 (2010) 
[8] F.N. Loreti, Y. Z. Qian, G. M. Fuller and A. B. Balantekin, 
Phys. Rev. D 52 6664 (1995) 
J. P. Kneller, arXiv: 1004. 1288 [hep-ph]. 
[10] Fogli, G.L., Lisi, E., Mirizzi, A., & Montanino, D., Phys. 
Rev. D 68 033005 (2003) 

J. P. Kneller and G. C. McLaughlin, Phys. Rev. D 73 
056003 (2006) 

J. P. Kneller and G. C. McLaughlin, Phys. Rev. D 80 
053002 (2009) 

[13] Kmetic, M.A., & Meath, W.J., Phys. Lett. A 108 340 (1985) 
[14] L. Wolfenstein, Phys. Rev. D 17 2369 (1978) 
[15] S. P. Mikheev & A. I. Smirnov, Nuovo Cimento C 9 17 
(1986) 



[4] 

[5] 
[6] 



[9] 



[11] 
[12] 



